Intra- molecular refrigeration in enzymes 
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We present a simple mechanism for intra-molecular refrigeration, where parts of a molecule are 
actively cooled below the environmental temperature. We discuss the potential role and applications 
of such a mechanism in biology, in particular in enzymatic reactions. 



Introduction. In the present paper, we would like to 
suggest the possibility that enzymes can cool their active 
sites below the environmental temperature and thereby 
increase their functionality. A similar cooling mechanism 
might be used by molecular machines to cool parts of the 
machine and to reduce the deteriorating effects of ther- 
mal noise. The very concept of intra-molecular refriger- 
ation was proposed by us in a recent paper [1]. Here, we 
discuss it in more detail and present, for the first time, a 
possible molecular mechanism. 

Let us first discuss why be believe that molecular cool- 
ing is important and what functionality it could have in 
biological systems. To start with, we note that for many 
species, including mammals and birds, the ability to cool 
parts of their body is essential to survive. Their body 
temperature is kept largely stable at a specific temper- 
ature even when they live in an environment with vary- 
ing temperatures. This is only possible since they have a 
built-in cooling system. Such a cooling mechanism works 
however at a large scale - essentially at the scale of the 
whole multicellular organism and generally involves the 
interaction of many organs. The type of cooling we are in- 
terested here is, however, completely different and could, 
in principle, exist in all three domains of life - bacteria, 
archaea and eukarya. We suggest the possibility that 
cooling below environmental temperature may occur at 
molecular level. 

If cooling could be achieved at molecular level, this 
would have obvious benefits. An example is an increased 
efficiency of catalysis: Many proteins act as catalyzers 
(enzymes) with very high specificity . They have active 
sites in the shape of cavities, which bind only molecules 
that fit precisely into the cavity, like a hand in a glove 
[see Fig. [ija)]. When a reactant molecule (substrate) 
binds to the active site, its activation energy to react 
with some other molecule is lowered, and the reaction is 
thus sped-up due to the presence of the enzyme. After 
the reaction has taken place, the reactant molecule (prod- 
uct) leaves the cavity. The temperature dependence of 
the reaction rate in the presence of an enzyme typically 
looks like the curve displayed in Fig. [ijb). At moder- 
ate temperatures, the reaction rate increases with the 
temperature through an increased thermal energy of the 
reactant molecules. At higher temperature, however, vi- 



brations in the protein may lead to a deformation of the 
cavity, which leads in turn to a reduced efficiency to bind 
reactant molecules (and eventually to de-naturation) and 
hence to a decreased rate of reaction. From this descrip- 
tion [2 it is clear that, if an enzyme would succeed in 
cooling its active site below the environmental tempera- 
ture by no-matter-what mechanism, then the enzymatic 
reaction will become more resistant and functionality will 
be increased. Cooling could, in principle, be achieved by 
different methods. The model we describe here is based 
on conformational changes, which is a common feature 
in molecular dynamics. 
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Figure 1: (a) Molecular model of catalysis: Binding of a re- 
actant molecule (blue) to the active site of the enzyme, (b) 
Temperature dependence of a chemical reaction which is cat- 
alyzed by an enzyme. 

Molecular mechanism. Let us now describe the 
molecular model in a little more detail. As we have men- 
tioned, the basic mechanism that is responsible for catal- 
ysis is the binding of the substrate to the active site of the 
enzyme, which has the effect of lowering the activation 
energy for the chemical reaction. 

The shape of the cavity at the active site is defined 
by a specific arrangement of functional groups, which 
matches with the molecular structure of the substrate. 
When the ambient temperature increases, there will be 
more frequent (and higher energetic) collisions of reac- 
tant molecules with the enzyme which increases the re- 
action rate. However, when the temperature becomes 
too high, it will lead to increased thermal vibrations of 
the functional groups near the active site as indicated in 
Fig. [2|a). As a consequence, the substrate will no longer 



be able to bind to the enzyme and the reaction stops. 

The simple refrigeration mechanism we propose is 
based on a conformational change of the enzyme, which 
contracts the cavity at the active site and thereby stabi- 
lizes its vibrating parts. The reason is that, if the cavity 
contracts, its walls will interact with each other, mak- 
ing the whole cavity "stiffer". Suppose the cavity stays 
closed until it reaches thermal equilibrium. At equilib- 
rium, due to their mutual interaction, the amplitude of 
the movement of the walls relative to each other is re- 
duced as compared to what it was when the cavity was 
open, see Fig. [lib). Indeed, we can see this in a differ- 
ent way: if the amplitude of the oscillation in the stiffer 
configuration would be the same as in the open configu- 
ration, it would cost more energy, i.e. above the average 
energy corresponding to the ambient temperature. Sup- 
pose now the cavity re-opens at a speed fast relative to 
the thermalization rate, but slow compared to the period 
of vibration. In this case, the opening does not itself am- 
plify the amplitude of vibration. Hence the cavity ends 
up in the open position with reduced vibrations, i.e. at 
a lower temperature, as is indicated in Fig. [2|c). The 
cavity is now cold and ready to accept new react ants. 

As in any refrigerator, in order for cooling to work, 
there must be a supply of free energy. We suggest as 
a possible mechanism that the conformational change is 
induced by an allosteric process, in which the effector 
molecule supplies the free energy, as illustrated in Fig. |2j 
However, we emphasize that the essential part of the pro- 
posal is the refrigeration via conformational change. How 
exactly the conformational change is produced, and how 
the free energy is delivered, is irrelevant. 

Experiment. An experimental demonstration of 
molecular cooling could be straightforward. Suppose one 
finds a catalyzing enzyme with an active site as described. 
Several measurement rounds are then performed, observ- 
ing the enzyme activity (reaction rate) as a function of 
temperature. In each round, the same concentrations of 
enzyme and substrate are used, but the concentration 
of the effector molecules is varied. The predicted tem- 
perature dependence of the chemical reaction rate would 
then be as described in Fig. [sj^a) (red curve) . The enzyme 
will retain its function up to higher temperatures due to 
the molecular cooling effect when effector molecules are 
present, compared to the situation when there are no ef- 
fectors. It is also possible that for lower temperatures 
the cooling may have a detrimental effect on the enzyme 
efficiency, since during the time when the active site is 
closed, it cannot accept any substrate, see Fig. |3|b). 

Physical model. Although we talk about cooling at 
a molecular level, the refrigeration mechanism is essen- 
tially classical, that is, not quantum mechanical. For our 
purposes, the different parts of the enzyme constituting 
the active site can be viewed as mechanical systems of 
different masses, which can oscillate around their equi- 



2 



(a) 




Figure 2: Illustration of refrigeration via allosteric conforma- 
tional change, (a) At very high temperatures, deformation 
and vibrations of the cavity at the active site will prevent re- 
actant molecules to bind to the enzyme, (b) Cooling of the 
active site can be achieved through a conformational change 
that contracts the cavity temporarily. Since in the "closed" 
position the cavity is stiffer, at thermal equilibrium the ampli- 
tude of vibrations is smaller than in the "open" position. Free 
energy needs to be supplied, e.g. by ATP or a similar effector 
molecule, docking to another active (so-called allosteric) site, 
(c) Upon reopening of the cavity, the amplitude of vibration 
remains reduced, until re-thermalisation. 



librium configuration. At increasing temperatures, the 
different parts of the system start vibrating with increas- 
ing amplitudes, which diminishes the functionality of the 
enzyme. 

The active site of the enzyme can vibrate in many 
complicated ways. Technically, there are many vibra- 
tion modes. For simplicity, here we analyze the case of 
just two oscillators, which stresses the essential features 
of the model. The two oscillating parts could represent 
the two "lips" of the cavity in Fig. [2|a). The vibrations 
indicated in the figure then correspond to excitations of 
the oscillators around their respective equilibrium posi- 
tions. Contracting the cavity brings the parts together; 
when the two lips are closer, additional forces will act on 
them. Whatever the sign of these forces (attractive or 
repulsive) will be, the net effect is to make the entire sys- 
tem tighter and increase the frequency of the vibrations. 
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Figure 3: Proposal of an experiment, (a) Predicted tempera- 
ture dependence of the rate of an enzymatic reaction, without 
and with (red) supply of effector molecules, (b) For lower tem- 
peratures, the cooling may have a detrimental effect on the 
enzyme efficiency. 



This process can be modeled by an extra interaction 
potential which is switched on in the contracted config- 
uration. While the exact nature of these interactions is 
not essential, its main effect is to increase the frequency 
of the oscillators which will be energetically more difficult 
to excite. 

The Hamilton function for this system is given by 



H{t) = 
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(1) 

where Xi is the displacement of mass from its equi- 
librium position, Pi the corresponding canonical momen- 
tum, and Ki the spring constant associated with the ith 
oscillator. The last term describes the variable interac- 
tion potential, which is switched on in the contracted con- 
figuration, accounted for by the variable spring constant 
K(t). Since we are only interested in the system near 
the equilibrium, we can choose quadratic oscillator po- 
tentials for all the involved potentials. For simplicity, we 
shall assume equal masses and spring constants for both 
oscillators, i.e. mi = m2 = m and ki = K2 = ^- It is con- 
venient to introduce new variables, Xcm = (-^1 +-^2)/2, 
^cm = Pi-\-P2 for the center of mass, and x = Xi — X2 , 
p = (P2 — Pi)/2 for the relative displacement coordinates 
and momenta. The transformed Hamilton function then 
reads 
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with /^cm = 2/^, Krei = M = 2m, /i = m/2. The mo- 
tions of the centre of mass and the relative displacement 
coordinates are decoupled and describe harmonic oscilla- 
tions. While the first does not lead to a deformation of 
the cavity, the latter does. We will thus concentrate on 
the motion of the relative displacements in the following, 
described by 
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Krel + K{t) 2 
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(3) 



When the cavity is open, ^{t) :^ 0, and, due the influ- 
ence of the environment, the relative mode of oscillation 
will be thermally excited. Upon contracting the cavity, 
the coupling function n{t) will effectively be switched on 
and tighten the potential. This corresponds to an in- 
crease of the spring constant, which makes this oscillator 
stiffer. To understand how this can lead to cooling, it 
is convenient to use a quantum mechanical description 
of the oscillator, where the energies are quantized. We 
would like to emphasize however that we consider a situ- 
ation that is essentially classical, that is, we are far away 
from the quantum mechanical ground state. Using the 
quantum mechanical language is here simply a technical 
convenience; we could also have solved the system clas- 
sically. 

The Hamilton operator corresponding to ([3| can be 
written in the form 



H(t) = nw{t) a^a 



with a time-dependent frequency 



and ladder operators 

a = at = 
t t 



(4) 

(5) 

(6) 
(7) 



For a fixed frequency, uo{t) = uo = const, the possi- 
ble energies of this oscillator are given by the formula 
En = hio{n + 1/2), n = 1,2,.... In a thermal environ- 
ment at temperature T, the different energies En will 
be found with probabilities Pn = ex.-p{—nhuj/kBT), 
with kB and h denoting the Boltzmann and Planck con- 
stant, resp., and Z a normalization factor (called the par- 
tition function) which ensures Xl^Pn = 1- 

In the open and contracted configuration, the oscilla- 
tor will have two different frequencies, loq and coi > cjq, 
respectively. Suppose that the system starts at time to 
in the open configuration, with a spacing of the allowed 
energies AE^ = hcoo^ and with a thermal-state distribu- 
tion j9n(^o) = ex.-p{—nhLJo/kBT). Upon contraction, 
the spacing between the allowed energies will increase 
to AE^^ = hijJi > hcjQ (see Fig. [4| . The thermalization 
of the oscillator is a process with a characteristic time 
scale that depends on the interaction strength between 
the system and the environment. Ideally we assume 
that, after contraction, the system stays long enough 
in this configuration, such that it has time to thermal- 
ize. The resulting thermal distribution at time ti is then 
Pn{ti) = Z~^ ex.-p{—nhuJi/kBT). Note that, since the 
spacing between the energy levels has increased, a larger 
fraction of the population will move to lower energies, as 
indicated in Fig. |4j 
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Figure 4: Refrigeration of the molecular vibrations through a 
conformational change (see text). 



Upon re-opening, the system will move back to the 
open configuration (with frequency cjq). Let the time re- 
quired for opening be Topen- Whether or not cooling takes 
place depends on two time scales. The first is the ther- 
malization time rtherm and the second the period of oscil- 
lation Tosc = 27rc<;~^. The ideal regime is when the tran- 
sition is fast compared to thermalization but slow com- 
pared to the oscillation period, i.e. Tosc <^ ''"open <^ '?"therm- 
In that case the system will cool down. Indeed, in this 
regime the opening is adiabatic and the relative popula- 
tion of the different energy levels does not change during 
the transition. In this ideal case, right after returning 



open? 



to the open configuration, i.e. at time ^2=^1+ r( 
the population of the energy levels will still be the same 
as at time ti, that is Pn{t2) = ex.-p{—nhu;i/kBT) 
= ex.-p{—nhuJo/kBT'^)' With respect to the new en- 
ergy spectrum (spacing AE^ = hujo)^ this corresponds 
however to a thermal state with reduced temperature 
= ^T. After such a cooling cycle, the temperature 
of the active site of the enzyme has thus been reduced by 
the amount AT/T = 1 — uq/uji. We note however that, 
as shown later, the cooling effect takes place, albeit with 
reduced efficiency, also outside the ideal regime described 
above. 

In the remainder of this paper, we illustrate the es- 
sential dynamics of the refrigeration process for a time- 
dependent oscillator frequency uj{t). For the present pur- 
pose, it suffices to describe the thermalizing effect of the 
environment by a master equation of the form 

^^p{t) = -'-[H{t),p{t)]+L,p{t) (8) 
with H{t) as in Q and with a dissipative term 



Ltp(t) = — — (i^ + 1) {a^ap + pa^a — 2apa^} 
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V {aaV + pcio) — 2a^ pa] . 



(9) 



Here p(t) is the state (density matrix) of the system at 
time t. The parameter 7 is the thermal relaxation rate 
that describes how fast the system would approach, for 
any fixed configuration {uo = const), its thermal equilib- 
rium state 



Ptherm(^) p{t ^ Oo)|c, 



exp 



a 



and V is the mean number of thermal excitations in that 
equilibrium state, v = u{uj) = l/[ex.p{huj/kBT) — 1]. For 
a time-dependent oscillator frequency, u = uj{t), this 
"instantaneous" thermal state thus changes with time, 
as does u = iy[u;{t)] . Unless the motion of the oscil- 
lator is very slow compared to the relaxation time, it 
will be driven out of equilibrium, which gives rise to the 
described refrigeration effect. In the regime where the os- 
cillation period is short compared to all other timescales 
(adiabatic regime), Tosc <^ ''open, ''"therm, the solution can 
be written in the form 3 
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where ri{t) is a solution of the differential equation 
d 



dt 



7jit) = --/7j{t)+-/Mu;it)] + l} 



(11) 



(12) 



which can be readily integrated to give 



^(t) = e-^^'-'Mh) +1 fds -^—^ . (13) 

Jtx 1 — e 

The function ri(t) has a simple interpretation, it is equal 
to the mean number of excitations (plus unity) in the 
oscillator at time t, i.e. r]{t) = {a) a) + 1. The state (11) 



describes a "quenched" Boltzmann distribution, corre- 
sponding to a quasi-thermal state, albeit with a tem- 
perature T{t) reduced below the temperature T of the 
environment. 

Tit) _ log Ht)/{i.it) + 1)] _ log Ht)/ii.{t) + 1)] 



T log[l-l/,?(t)] log[(ata)/((ata) + l)] ■ 

(14) 

Let us now have a closer look at the essential part 
of the process, that is, when the molecule moves from 
the closed configuration (after thermalization) at time ti 
back to the open configuration, which it reaches at time 
t2 = ti-\- Topen- In Fig. [sj we plot the temperature of the 
oscillator - representing the active site - as a function of 
time. As an example, we choose the specific time profile 
uj{t) =cji + (cJo-^i)(sin(f :^)), for = ti < t < ^2 = 
Topen and uj{t) = ujQ for t > Topen, but othcr profiles give 



similar results. The function r]{t) in (13), like u;{t), can 



(10) 



be expressed in terms of the dimensionless quantities ^ , 

A crucial issue is, of course, to determine the ranges of 
values of the parameters (opening time, oscillation fre- 
quencies and thermalization time) relevant to our sys- 
tem. This is however not so straightforward, since these 
parameters vary over considerably large ranges [4] . They 
depend on the type of molecule we consider, on its size, 
on the size of the active site, on which parts of the ac- 
tive site contract (i.e. the whole site or only subparts 
involving some of the functional groups), and on the de- 
gree of interaction with the environment. Furthermore, 
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such parameters are typically not all available for the 
same molecule, but some parameters are known for some 
molecules, other parameters for other molecules. 

For the figure, we have chosen the following values: 
= 0.032, ^ = 2, 7ropen = 1- At room tempera- 
ture, T = SOOiC, this corresponds to molecular oscillation 
times Tosc = 27r/a;o = 5ps, r^g^ = 2i\ jujx = 2.5ps, which is 
a typical time scale of elastic vibrations. Possible values 
for the opening time Topen, which is the time for (part of) 
a conformational transition, may range from fractions of 
nanoseconds up to microseconds and beyond, depending 
on the size of the protein and the specific type of tran- 
sition. Similarly, the thermalization time Ttherm = 
may vary significantly, depending on coupling of the os- 
cillators to the environment. For fast configurational pro- 
cesses, an exemplary choice would be Topen = lOOps as 
the time scale for re-opening the cavity and = lOOps 
as the thermalization time for the vibrations. 

As we have mentioned, these values will depend on 
the specific molecular realization, and all of these num- 
bers are subject to considerable variation [5 . But it 
seems reasonable to assume that the time scales of the 
conformational transition and of the cavity vibrations 
are separated by at least one order of magnitude, i.e. 
'^osc, '^osc ^ '^open (which is the so-called adiabatic regime, 
with uj/uj <C 1). The timescale for Topen may also be 
much larger; likewise, the timescale for thermalization 
may be different. Whether or not the refrigeration ef- 
fect can be observed depends on the timescale of the 
thermalization relative to the opening time, i.e. on 
'Topen /Ttherm = T^open- As loug as the thcrmalizatiou 
time is larger or at least comparable to the opening time, 
i.e. 7Topen ^ 1, the effect can clearly be observed, as 
demonstrated in Fig. [5] One can seen how the effective 
temperature of the cavity T{t) drops during the opening 
transition to about 65% of the environmental tempera- 
ture T. After the cavity has reached the open position, 
its temperature relaxes to the environmental tempera- 
ture. It should however be noted that a useful effect of 
refrigeration may persist up to several multiples of Topen- 
A refrigeration by about 1 degree Kelvin, for example, 
means a relative value T{t)/T = 299/300 = 0.997 which 
is reached only after about 6 x Topen- Of course, the effect 
is more persistent for longer thermalization times. 

Conclusions. The mechanism for intra- molecular re- 
frigeration that we have discussed in this paper is simple. 
It is based on conformational changes that induce time- 
dependent forces between different parts of a molecule, 
reducing their relative vibrations. It works as long as the 
configurational changes, which may involve only small 
parts of a protein, happen faster or at least on a compa- 
rable timescale as the damping of the relevant vibrations. 
The changes in local temperature, which can be gener- 
ated by this mechanism, are, even for moderate variations 
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Figure 5: Refrigeration of the active site as a function of 
time, (a) Temperature T{t)/T of the active site (blue) and 
frequency of vibrations a;(t)/a;i (red), (b) Average number of 
quanta {a^a){t) in the vibration mode. [Parameters: ooi/ooo — 
2, 7ropen = 1, ^ = 0.032. At T = 300i^, this corresponds 
to Tosc = 2n /ujQ = 5ps as the time scale of elastic vibrations.] 

of the binding potential, quite significant: It should be 
remembered that, for many biological processes, a refrig- 
eration by only a few degree Kelvin can make all the 
difference. 

Many variations of this scheme are conceivable. For 
example, it is possible that cooling due to conformational 
variations may take place at other sites, near the active 
site, but without changing the latter's shape. Through 
heat transfer to the cooling centers, effective refrigeration 
of the active site would be possible without interrupting 
the catalyzer process. 

The improvement of enzyme functionality is but one 
example of the potential applications of intra-molecular 
refrigeration, but there should be many others. Molecu- 
lar refrigeration might also play a role e.g. in molecular 
machines, to stabilize their operation against the accu- 
mulation of noise and errors, or in more complex pro- 
cesses such as those involved in protein synthesis. 

Last, but not least, the possibility of maintaining tem- 
perature gradients on a molecular scale is by itself of 
considerable interest for the investigation of thermody- 
namical processes at the bio-molecular level. 
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